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Abstract 



A major mystery of glass- forming liquids is the non-Arrhenius 
temperature-dependence of the average relaxation time. This paper briefly 
reviews the classical phenomenological models for this phenomenon - the free- 
volume model and the entropy model - and critiques against these models. We 
then discuss a recent model [Dyre, Olsen, and Christensen, Phys. Rev. B 53, 
2171 (1996)] according to which the activation energy for the average relax- 
ation time is determined by the work done in shoving aside the surrounding 
liquid to create space needed for a flow event. In this model the non-Arrhenius 
temperature-dependence is a consequence of the fact that the instantaneous 
(infinite-frequency) shear modulus increases upon cooling. 
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I. INTRODUCTION 



Apparently all supercooled liquids are able to form glasses fIHT2fl. The glass transition 
takes place when the viscosity of the supercooled liquid upon cooling becomes so large 
that molecular motion is arrested on the time-scale of the experiment. The fascination 
of this phenomenon lies in the fact that chemically very different liquids - involving ionic 
interactions, van der Waals forces, hydrogen bonds, covalent bonds, or even metallic bonds 
- exhibit a number of common properties when cooled to become highly viscous P,|S|,|T|,|T2"|| . 
Of particular interest here is the temperature-dependence of the average relaxation time, 
t. This quantity may be determined, e.g., as the inverse dielectric, mechanical or specific 
heat loss peak frequency. Alternatively, it may be calculated from the viscosity rj and the 
infinite-frequency shear modulus by means of Maxwell's expression 

r = £- ■ (1) 

These definitions do not give exactly identical r's, but the difference is insignificant for the 
present purposes. It is widely believed [p],^|,^3|- |2^j31| that different r's are roughly identical 



because they basically measure the rate of "flow events" : Most molecular motion in a highly 
viscous liquid is purely vibrational around a potential energy minimum. Only seldom does 
real motion take place. This happens in the form of a sudden rearrangement of molecules, 
a process which is unlikely because of the large potential energy barrier to be overcome 
|l| Jl5|j22]| . Kauzmann referred to flow events as a "jumps of molecular units of flow between 



different positions of equilibrium in the liquid's quasicrystalline lattice" [Q. The molecules 
involved in a flow event define a "relaxing unit" [0, "cooperatively rearranging subsystem" 



or "cooperatively rearranging region" |y|], "quasi-independent unit" | 17| , "thermokinetic 



structure" p0| , "molecular domain" p2fl , or "dynamically correlated domain" j23 



As the glass transition is approached, the average relaxation time becomes longer 
and longer. For typical cooling rates r is of order 10 3 s at T g . From a general physi- 
cal/chemical point of view, the temperature-dependence of r is anomalous in the following 
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sense. In only very few liquids is r Arrhenius (examples are SiC>2, GeC>2, BeF2 or albite 
(NaAlSiaOg) 0). Predominantly, r is non-Arrhenius by exhibiting an apparent activation 
energy [d In r / d{kBT)~ 1 \ that increases as the temperature decreases . A measure of the 
departure from non-Arrhenius behavior is the fragility m, defined as the apparent activation 
energy at T = T g in units of ksT g ln(10) For a simple Arrhenius liquid m is about 16; 
for most viscous liquids m is between 50 and 150. Liquids with large fragility are termed 
"fragile", liquids with fragility not far above 16 are termed "strong" |2q| . There is a gen- 



eral tendency that fragile liquids have broader distributions of relaxation times than strong 



liquids [p7| , Pq| . This rule, however, is not without exceptions p9| , |30 . 

In the discussion below we will not distinguish between the temperature-dependence 
of average relaxation time and of viscosity, because these two quantities are roughly pro- 
portional (in Eq. ([]]) the temperature-dependence of is insignificant). We identify the 
temperature-dependent activation energy from the expression |]31|j3^| 



T 



Although AE(T) is different from the apparent activation energy, experiment imply that 
AE(T) also increases as the temperature decreases. 

It is not at all obvious that a general explanation for the non-Arrhenius behavior of chem- 
ically quite different viscous liquids exists, but it seems to be a reasonable first hypothesis. 
This paper discusses models for the non-Arrhenius average relaxation time, models that are 
phenomenological in the sense that r(T) is determined by some macroscopic property of 
the liquid. The most famous phenomenological models for the non-Arrhenius r's are the 
free-volume model of Cohen, Turnbull and Grest |I^ , |33| - p5|1 and the entropy model of Gibbs, 



DiMarzio and Adam [ 13 , 36 1 . These models and critiques against them are briefly reviewed 



below (see also Johari's review of phenomenological models Ref. [ffl). We then discuss a 



recently proposed model |37],|38|], according to which the activation energy of a flow event 
mainly originates in the work done in shoving aside the surrounding liquid to create enough 
space for the molecules to rearrange. 



II. EARLY PHENOMENOLOGICAL MODELS 



The importance of volume was stressed long time ago by Eyring and coworkers, who 
suggested that the viscosity of a liquid is lower the greater the number of holes present |39| . 



Defining the free volume per molecule Vf as the average volume per molecule in the liquid 
minus the volume of the molecule itself, Doolittle |40| in 1951 found that the viscosity of a 



number of simple hydrocarbon liquids may be fitted by the expression 



rj = r] exp ( — j . (3) 
\ v fJ 



In 1959 this expression was derived by Cohen and Turnbull arguing as follows [Q. The 
molecules are modelled as hard spheres. A molecule is mostly confined to a cage bounded 
by its immediate neighbors. Occasionally, there is a fluctuation in density which opens 
up a hole within the cage. Molecular transport occurs only when a void having a volume 
greater than some critical value v* forms. The total free volume may be distributed in 
various ways between the cages. The average relaxation time is essentially the inverse of the 
probability P that redistribution of free volume by chance creates a void of greater volume 
than v*. Turnbull and Cohen calculated this probability by standard statistical mechanical 
arguments [|33|]. Their result is P cx exp(— C/vf), leading to Eq. (|||) via Eq. (fj). A basic 



assumption in the free-volume model is that no energy is required for redistribution of free 
volume. When the model is applied to real liquids, Cohen and Turnbull defined the free 
volume as that part of the excess volume that may be redistributed with no increase in 
energy 
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In the free-volume model the temperature-dependence of average relaxation time comes 
from the fact that the free volume decreases with decreasing temperature. If the free volume 
is taken to be a linearly decreasing function of temperature, one arrives at the famous Vogel- 
Fulcher-Tammann (VFT) expression, 

r = r exp (——] . (4) 



vT-T 0/ 

Here T is the temperature at which there is no free volume 



It is noteworthy that in the approach of Cohen and Turnbull, the concept of free volume 
has a meaning different from that of Doolittle ||40|| . In Doolittle's definition, the molecular 



volume is obtained by extrapolating the liquid volume to zero temperature and consequently 
the free volume is zero only at zero temperature. Thus, while the experimentally motivated 
Eq. (|3p was the starting point of Cohen and Turnbull, their theory represents a quite different 
way of thinking than that of Doolittle. 

What critiques may be raised against the free-volume model? The derivation of Eq. ([3]) 
may be questioned because of the primitive way in which the entropy of the free volume 



is taken into account. Moreover, despite several attempts |[14| , |34| the very concept of free 
volume in Cohen and Turnbull's sense seems to be ill defined operationally for general 
liquids. When it comes to a comparison to experiment, the VFT-equation Eq. (^) often 
gives a good fit to data |7]||. However, the fit is seldom perfect; in particular, there are 
systematic deviations close to T g , where the average relaxation time is apparently always 



less temperature-dependent than predicted by Eq. (f|) p| j3^j4T| -P4| . Since Eq. ([J) is derived 
by combining the free- volume model with the ad hoc postulate that the free volume depends 
linearly on temperature, a more direct test of the model may be performed by applying 
pressure to the liquid. The model predicts that the average relaxation time is solely a 
function of density. Indeed, r does increase dramatically at high pressures, but quantitatively 
the free- volume model is not confirmed f3[. A further test of the free- volume model is based 
on the fact that the glass transition is similar to a second-order phase transition in the sense 
of Ehrenfest (with continuity of first derivatives of the free energy and discontinuity of the 
second derivatives). For the pressure- dependence of the transition temperature, the claim 
that the average relaxation time is controlled by volume translates into the requirement that 



the glass transition takes place at constant volume. This implies E5] 



dTg AK 

dp Aa ' 1 ' 

where An is the difference between the isothermal bulk compressibility of liquid and glass 
and Aa the same difference for the isobaric thermal expansion coefficient. Equation (|5|) is 



seldom fulfilled fl46|j47fl. 



In the free-volume model the glass transition occurs at a fixed volume. In the theory 



of Gibbs and DiMarzio from 1958 |36| the variable controlling the average relaxation time 
is the configurational entropy. Evaluating the partition function for a lattice model of 
linear polymeric chains in a mean-field approximation, Gibbs and DiMarzio found that 
there is a second-order phase transition at a finite temperature Tk to a low-temperature 
state of zero configurational entropy. This state is a "ground state" of amorphous packing. 
Furthermore, Gibbs and DiMarzio argued that in the neighborhood of Tk the energy barrier 
restricting transitions between different molecular configurations is very high, because "the 
few states that could conceivably occur close to T K are widely separated in phase space, so 
proceeding from one to another involves a considerable change in the topology of molecular 
entanglements" . In this picture, the very equilibrium properties of a supercooled liquid give 
rise to kinetic sluggishness which prevents the equilibrium second order phase transition 
from being reached in finite time. 

These ideas were quantified in 1965 by Adam and Gibbs They argued that the size 
of cooperatively rearranging regions, defined as "the smallest regions that can undergo a 
transition to a new configuration without a requisite simultaneous configurational change 
on and outside its boundary", diverges as the configurational entropy goes to zero. The 
region size is estimated by requiring that at least two different configurational states should 
reside in a region, leading to a size inversely proportional to the configurational entropy, S c . 
If the energy barrier to be overcome is assumed to be proportional to region size, Adam and 
Gibbs arrived at the following expression for the temperature-dependence of the average 
relaxation time, 

t = r exp f — . (6) 



.S C T 

Close to Tk the denominator S C T may be expanded to first order in T — Tk whereby Eq. 
([]) becomes the VFT-expression Eq. (|) with 

T = T K . (7) 
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The entropy model resolves the Kauzmann paradox of a negative configurational entropy 
below Tk without just avoiding it (as Kauzmann did himself |1J by suggesting that crys- 
tallization sets in before Tr- is reached). In many respects the model is in good agreement 
with experiment. Thus, Eq. (|7|) is often obeyed P JT^r47|J^] , which is quite remarkable, 
given the dynamic definition of T and the quite different purely thermodynamic definition 
of Tk- In particular, systems with only small excess specific heat (implying less pronounced 
Kauzmann paradoxes and TVs close to zero), generally tend to be "strong", i.e., have VFT 
T 's close to zero [|^. 

The entropy model undoubtedly presents a beautiful scenario. Still, it may be critiqued 
both in regard to its relation to experiment and in regard to its inner consistency. Experimen- 
tally, there is no proof that a second order phase transition to a state of zero configurational 
entropy is underlying the laboratory glass transition. In many cases a simple two-level sys- 
tem fit excess entropy data well []50[ ; more generally, the data may be fitted with a model 
with only few energy levels fl?|| . Also, as mentioned above, the VFT-expression Eq. (|J) 
fails close to T g , where data are usually less temperature-dependent. Finally, it should be 
mentioned that the identification of excess entropy with configurational entropy rests on an 
assumption that the glass has the same "fast" contribution to the entropy as the crystal 
at the same temperature. As pointed out by Goldstein |5l]| , this assumption is not always 



realistic because the glass may have significant contributions to the "fast" specific heat from 
anharmonic vibrations and secondary relaxations not present in the crystal. 

In regard to the inner consistency of the entropy model, we first note that the mean- 
field solution of the lattice polymer model of Gibbs and DiMarzio is incorrect and that, in 



fact, the model has a positive configurational entropy at all positive temperatures ||52| , |53 



Ignoring this objection and accepting the general idea of a phase transition to a state of zero 



configurational entropy, one may reasonably ask What is the nature of the amorphous 
ground state, the "ideal glass state"? Since this state is unique a simple description of it 
would be expected; however none has been proposed. The argument of Adam and Gibbs 
is also not compelling. They assumed ad hoc that the energy barrier to be overcome is 



proportional to the size of the cooperatively rearranging region. Though this may seem 
reasonable, it does not have to be correct. More generally, approaching a zero-entropy state 
does not in itself imply a diverging relaxation time. There is no compelling link between 
dynamics and thermodynamics: In a master equation description of the dynamics many 
different possible forms of transition rates - leading to quite different relaxation behaviors - 
are consistent with the same statistical mechanics. 

Instead of focussing on volume or entropy as the variable controlling the relaxation of 
viscous liquids, potential energy may be the relevant variable, as first suggested by Gold- 
stein ||15|| . A number of authors have taken this approach P, p^JT7| , [2T| --p^, pl| , ^ , ^5|1 . In the 



simplest energy controlled models the transition state of a region is taken to be temperature- 
independent with potential energy E |16| , |17| , pT| , leading to the following expression for the 
activation energy in terms of the average potential energy of one region E(T) is arrived at 



0: 



AE(T) = E -E(T). (8) 

Since the average potential energy decreases with decreasing temperature, the activation 
energy increases. Qualitatively, this is what is seen in experiment. However, in order to 
fit data relatively large regions are needed, implying much broader relaxation time distri- 



butions than observed |31 |. Therefore, the simple picture does not work and more involved 



approaches need to be taken |55 . 



III. SHOVING MODEL 

The models above discussed all assume that relaxation depends only on the state of the 
region involved. A completely different approach may be taken, where the relaxation rate 
depends only on properties of the surrounding liquid: As starting point we take the fact that 
molecular interactions are strongly anharmonic, i.e., with strong short-ranged repulsions and 
weak long-ranged attractions; as shown by Chandler, Weeks and Andersen |5B| , this "van 
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der Waals" picture explains a number of phenomena in liquids. Next, as in the free- volume 
model we assume that space is needed for molecules in a viscous liquid to rearrange. The idea 
is that, because of harsh intermolecular repulsion, rearrangement at constant region volume 
is excessively costly, so it is much easier for the molecules to spend energy on shoving aside 
the surrounding liquid. If the rearranging molecules constitute a sphere which changes its 
radius by Ar, the energy cost for expanding is A(Ar) 2 (because the surrounding liquid may 
be regarded as an elastic solid on the short time scale of a flow event). The energy barrier 
to be overcome inside the sphere is some function /(Ar), that varies strongly with Ar. 
Minimizing the total energy cost leads to 2AAr + /'(Ar) = 0. If the ratio between the 
"shoving" work and the "inner" barrier to be overcome is denoted by A, we find 

A(Ar) 2 1 din/ 

/(Ar) 2dlnAr ' 1 ' 

Because of the strong repulsions one expects this logarithmic derivative to be numerically 
much larger than one. Thus, the shoving work gives the dominant contribution to the energy 
barrier. For simplicity we ignore the "inner" contribution to the activation energy. 

To calculate the "shoving" work we use the fact that during the flow event the surround- 
ing liquid behaves as an elastic isotropic solid with bulk modulus and shear modulus 
Goo. These elastic constants are known to be much more temperature-dependent in viscous 
liquids than in simple liquids or solids (crystals or glasses). Both and Goo increases 
as the temperature decreases. The work done on the surroundings depends linearly on 
these constants, thus leading to the observed increase in activation energy with decreasing 
temperature. 

Actually, it is only the shear modulus that is important. To show this we refer to the 
theory of elasticity of isotropic media ||57|| , assuming that the activation volume is relatively 
small: If V is the volume of the cooperatively rearranging region it is assumed that AV << 
V. We identify the activation energy with the elastic energy stored in the surroundings 
when the volume of the region has expanded to V + AV. Remember that elasticity theory 



| 57j concerns the relation between the stress tensor and the strain tensor Uij. The latter 



is defined by 

Uij = ^ (diUj + djUi) , (10) 

where di = d/dxi and Ui is the i'th component of the elastic displacement vector u. For an 
isotropic solid the bulk and shear moduli K and G are defined ]57f by 



Gij = KuuSij + 2G (u^ - T^SijUu^J . (11) 
The equation for static equilibrium is 

3^ = 0. (12) 

Substituting Eq. ( |I~0D into Eq. ([□]) and subsequently Eq. ( |TTD into Eq. (|12|) leads to 

(V+^G)V(V-u) + GV 2 u = 0. (13) 

For a purely radial displacement V x u = and thus, via the well-known vector identity 
V x (V x u) = V(V ■ u) — V 2 u, we have V 2 u = V(V • u). When this is substituted into 
Eq. ( |13D one finds 

V(V-u) = 0, (14) 

implying that V • u = C±, where C\ is a constant. The displacement (which is radial) is 
found by solving V • u = r~ 2 d r (r 2 u r ) = C*i, leading to u r = C*2 r ~ 2 + Cir/3. The latter term 
diverges as r — > oo and thus C\ = 0. In conclusion V ■ u = 0, i.e., there is no compression 
of the surroundings during a flow event. 

If the radius of the region before the expansion is R and the change of radius is AR, we 
have since AR << R 

R 2 

u r = AR — (r > R) . (15) 

The energy density of an elastic solid is |57j \Ku\ + G (v,ij — ^5ijUu^ . Since uu = the 
energy density is given by GuijUij = G{u 2 r + u 2 ^ + u1 e ) (all mixed terms like, e.g., u 2 ^ are 
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zero because the displacement is purely radial). When Eq. (|15D is used in the definition of 
the strain tensor in polar coordinates, we get for the energy density 6G(AR) 2 R 4 r~ 6 . Thus, 
the elastic energy is given by 



/ 6G(AR) 2 R 4 r- 6 Unr 2 )dr = 8nG(AR) 2 R. 
Jr 



(16) 



Substituting V = 47rR 3 /3 and AV = AnR 2 AR into Eq. flilf ) we find, introducing the 
"characteristic volume" 



2(AVy 



3 V 



(17) 



for the activation energy (with G = Goo(T)) 



AE(T) = G^T) V c . 



(18) 



For the average relaxation time we thus have [37 



t exp 



G^T) Vc 
k B T 



(19) 



Interestingly, extended mode-coupling theory leads to an expression resembling Eq. (|TD|) . 
except that G^ is replaced by the zero-frequency bulk modulus p8]| . The prediction of Eq. 
(|l9l) was checked against experiment on a number of organic liquids in Ref. Wn \, assuming 
that the characteristic volume is temperature-independent (strictly speaking, this assump- 
tion is inconsistent with Eq. (|9|), but for a strongly anharmonic potential the temperature- 
dependence of V c is negligible compared to that of Goo). I n Ref. [0 the following version 
of the well-known "Angell plot" p6[ was used: Instead of plotting the logarithm of the vis- 
cosity as function of T g /T, it was plotted as function of x = Goo(T)/T (normalized to one 
at Tg). The model predicts that a straight line should result. In Fig. la of Ref. |37[ results 
are shown based on measurements of the frequency-dependent shear modulus covering the 
frequency range lmHz-50kHz. Figure lb of Ref. [|37] shows the data of Lamb and coworkers 
from 1967 ||5!|, where G^ was obtained by an ultrasonic standing wave technique operating 
in the MHz region. Finally, Fig. lc of Ref. |37] shows data for two liquids, where G^ was 
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obtained by transverse Brillouin scattering. Overall, we found good agreement with the 
model prediction. 

It is generally believed that initial stages of glassy relaxation proceeds with a 
temperature-independent activation energy, which is a characteristic of the frozen struc- 
ture. According to the present model, however, there is a slight temperature-dependence of 
the activation energy for glassy relaxation, deriving from the fact that Goo is not completely 
temperature-independent in the glassy phase. The fact that Goo is measurable means that, 
if the model is valid, it is possible to monitor the activation energy for glassy relaxation 
directly. In fact, if one assumes Eq. fll9|) for the non-equilibrium structural relaxation time, 
the model gives definite predictions for the rate of relaxation of Goo itself. Note that, because 
Goo determines the relaxation rate of 

equilibrium liquid as well as of glass, this quantity is a direct measure of the Active 
temperature. 

IV. DISCUSSION 

Early phenomenological models link the non-Arrhenius average relaxation time of viscous 
liquids to configurational entropy, free volume, or potential energy. Here, an alternative 
approach to the non-Arrhenius problem was taken, linking t(T) to the high-frequency elastic 
shear modulus. The starting point of the shoving model is the fact that intermolecular 
forces are strongly anharmonic. Anharmonicity enters the model at three stages in the 
argumentation. First, the strong repulsions imply that it is very costly for molecules to 
rearrange at constant volume (a qualitative argument reminiscent of the free- volume model). 
Secondly, anharmonicity implies that the shoving work much exceeds the "inner" energy 
barrier (Eq. (|^), a quantitative argument). Finally, the fact that Goo depends on temperature 
is itself a consequence of anharmonicity. Note that, according to the model, one expects a 
liquid to be more fragile the more anharmonic it is. This, in fact, is what Angell conjectured 
arguing within the entropy model [|TD[] . 
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The shoving model basically involves three postulates: 1) The main contribution to the 
activation energy is elastic energy; 2) This elastic energy is located in the surroundings of 
the reorienting molecules; 3) The elastic energy is mainly shear energy. It is interesting to 
note that the model of F. Bueche from 1959 |Kj also focusses on the elastic properties of 
the surroundings: A particular molecule is regarded as surrounded by spherical shells of 
molecules, shells that are bound elastically to each other. Bueche's idea is now that if all 
concentric shells should vibrate outward in phase, the innermost shell would expand greatly, 
leaving the central molecule in a rather large hole so it could move to a new position. To 
calculate the probability of this happening Bueche made some further assumptions, leading 
to an expression that at high temperatures gives a simple Arrhenius expression but at low 
temperatures a VFT-expression. 

Returning to the shoving model, even if one basically accepts the above three postulates, 
there are a number of points potentially leading to deviations from Eq. fll9|): 1) Eq. flT9|) is 
based on a continuum approximation that may not be applicable on the molecular level; 2) 
The "inner" contribution to the activation energy has been ignored; 3) In real flow events 
spherical symmetry is probably violated to some degree, leading to some compression of the 
surroundings and thus a contribution to the activation energy proportional to K^; 4) In 
comparing Eq. (|i~9"D to experiment we have ignored any temperature-dependence of V c . 

Reference |37j gave a discussion of models related to the shoving model. To the best 
of the author's knowledge, the first to predict an expression equivalent to Eq. ([19]) was 



Nemilov [61] who - arguing quite differently - in 1968 arrived at this expression with our V c 



identified with the total region volume. At the present meeting Buchenau presented a model 



also leading to Eq. ( |T9|) for the average relaxation time |62j| . As emphasized by Buchenau 
both here and previously p3]| , models of this type link short time dynamics with long time 
dynamics. At first sight such a link may seem surprising but it makes sense 



because the transition itself is a very fast process |64|| . Finally, we note |k| that the 
present mechanism may possibly be applied also to explain the non-Arrhenius relaxation 
times of plastic crystals ||| and orientational glasses ]67| . 
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